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We propose and analyze neutral atom traps generated by vortices imprinted by magnetic field pulse 
sequences in type-II superconducting disks and rings. We compute the supercurrent distribution 
and magnetic field resulting from the vortices in the superconductor. Different patterns of vortices 
can be written by versatile loading field sequences. We discuss in detail procedures to generate 
quadrupole traps, self-sufficient traps and ring traps based on superconducting disks and rings. The 
ease of creating these traps and the low current noise in supercurrent carrying structures makes our 
approach attractive for designing atom chip interferometers and probes. 

PACS numbers: 37.10.Gh, 03.75.Be, 74.78.Na 



I. Introduction 

Recently superconducting atom chips have generated 
a lot of interest due to their attractive properties, such 
as the Meissner effect for type-1 superconductors and 
vortices for type-II superconductors [TH3]- Cold atoms 
have been trapped and manipulated near superconduct- 
ing surfaces [3H5] and Bose-Einstein condensation (BEC) 
on a superconducting atom chip has been reported re- 
cently |B]. Thermal and technical noise in proximity to 
superconducting surfaces have been shown both theoreti- 
cally and experimentally to be significantly reduced com- 
pared to conventional atom chips [7Ml4j. Superconduct- 
ing atom chips have the potential to coherently interface 
atomic and molecular quantum systems with quantum 
solid state devices [T5HT9] . Furthermore, ultracold atoms 
can be used to probe properties of superconductors, such 
as the dynamics of vortices in type-II superconductors 

EDI EH- 

Previously we reported on the design and realization 
of several types of magnetic traps involving linear type- 
II superconducting microstructures [22Tf24] . These traps 
are formed by the magnetic field carried by vortices im- 
printed in the superconductor to confine ultracold atoms. 
Loading of the vortices is accomplished by tailored exter- 
nal magnetic pulse sequences applied to the superconduc- 
tor in the mixed state. In a different approach vortex- 
based magnetic traps have been observed after cooling a 
superconducting disk through the critical temperature in 
the presence of an external magnetic field [21] . 

In this article, we extend our investigation of vortex 
based microtrap geometries to type-II superconducting 
disks and rings. To this end we have developed a numer- 
ical approach to calculate the vortex patterns imprinted 
after applying various external magnetic field pulses to 
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the superconducting structure. Using this numerical ap- 
proach we design quadrupole traps, self-sufficient traps 
and ring traps in the various superconducting geome- 
tries. We demonstrate that, by simple application of 
an additional bias field, a single quadrupole trap may 
be deformed into a ring. Considering the low current 
noise expected from the persistent supercurrents gener- 
ating these traps, and the low thermal noise expected 
from superconducting atom chips, this method may be an 
attractive way to generate low-noise on-chip atom inter- 
ferometry. Compared to other low noise traps generated 
by permanent magnets our system retains one desirable 
feature: The ability to control the trap parameters and 
geometry from shot to shot. 

The paper is organized as follows. Sec. [TTJdescribes the 
equations of motion for the sheet current density J in a 
type-II superconducting strip, disk and ring in an exter- 
nal magnetic field. To validate our numerical approach 
we compute the supercurrents in a superconducting strip, 
and compare them with analytical results [55]. For the 
two-dimensional structures considered in this paper no 
explicit analytical solution is known. In Sec. |III| we de- 
sign several types of traps for cold atoms by tailoring the 
magnetic field pulses which induce vortex patterns in the 
superconducting disk and ring. Sec. IV summarizes our 
results. 



II. Loading vortices 

Consider a type-II thin superconducting strip, disk and 
ring in a magnetic field perpendicular to their surface (xy 
plane) as shown in figjl] Here d denotes the thickness of 
the superconductors and a denotes the half width of the 
strip and the (outer) radius of the disk and ring. The 
strip is assumed to be infinitely long along y, so that we 
can restrict our simulation to the xz plane. Similarly, we 
can reduce our simulation to the rz plane for the disk and 
ring due to rotational symmetry. Throughout the paper 
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FIG. 1: (Color online) Schematic of the investigated type-II 
superconductor geometries. All geometries are assumed to be 
in the thin film limit with a thickness d. Shown are a strip of 
width 2a (left), a disk with radius a (center) and a ring with 
an outer radius of a and an inner radius of 0.4a (right). The 
magnetic field pulses are applied perpendicular to the surface. 



we will assume that the temperature of the superconduc- 
tors is below the critical value T c and that the thickness 
d is small compared to the other length scales of the 
system. When the external field B ext exceeds the lower 
critical field B c i, but lies below the upper critical value 
B C 2, magnetic flux starts to penetrate the superconduc- 
tor in the form of Abrikosov vortices [25]. To facilitate 
our numerical analysis we model these vortices together 
with the screening current as a macroscopic supercur- 
rcnt J(u,t). Our approach closely follows the method 
outlined in [26| I27j which makes use of the continuum 
approximation where the flux line spacing and the mag- 
netic field penetration depth are small compared to all 
other relevant length scales considered [35]. From the 
distribution of the supercurrent we then compute the to- 
tal magnetic field outside of the superconductor using the 
Biot-Savart law. We approximate the current density in 
the superconductor to be constant over its thickness d, 
which is valid in the thin-film approximation where d is 
smaller than the London penetration depth of the super- 
conductor. The sheet current density is then denned as 
J{x) = j(x)d, where j{x) is the local current density. 

The perpendicular field component B z (x) generated in 
the specimen plane z — by a constant B ext and by the 
sheet current J{x) along y for the strip can be written as 



B z {x) = B e 



+ ^[ a j(u)(^ i-)d«, (I) 

in Jo \ x — u x + uJ 



For cylindrically symmetric structures, where J(r) is cir- 
cling clockwise, we can write the perpendicular magnetic 
field component with ro = for disks and ro = 0.4a for 
rings as [2"o] : 
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where E(k) and K{k) are the first and second kind of 
the complete elliptic integrals with k = \/4ru/ (r + u) 
and /j,o is the vacuum permeability 2(i. 27] . In order to 
find an expression for the sheet current density J(x) we 
take into account the time-dependence of B ext and use 
Maxwell's equations, to arrive at the equation of motion 



for the sheet current density J(x, t) of the strip 
dB ext [ a dJ(u,t) Ar 
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dt Jo dt 
and J(r, t) for disks and rings [35] : 
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Here, p( J) — p c (J/ Jc)"^ 1 is the nonlinear resistivity with 
critical value p ci and n — Uo/k^T where Uo is the charac- 
teristic activation energy of the superconductor |29| . J c 
is the critical sheet current density [23111], defined as the 
maximum current density the superconductor can carry 
without a transition into the normal state. For the scope 
of this paper we have chosen the value of n = 19. This 
value gives a moderate response of the superconductor 
to applied fields. The exact value for n is material de- 
pendent and does not qualitatively alter the characteris- 
tic current distributions discussed for the experimentally 
applicable values. 



We approximate the integration in Eq(3] and [4] by a 
Riemann Sum. First, we discretise the width of the 
strip and disk (or ring) into N elements with the same 
size Ax — a/N or Ar = a/N(Ar = a — (ro)/N), and 
denote the middle point of each element by x n or r„, 
n = l,2...N. Then we rewrite Eqj3] in the form of a 
Riemann Sum in matrix form 



J{Xi,t) 



j—N 

dB ext v-^ nsr dJ{xj,t) 
2iTXi — h > M{xi,Xj) 



dt 



dt 



To solve the differential equation, we invert Eq(6]to find 
an equation tor — ~gt~ 



j=N 



= J2 M-^x^xj) 



J(xj,t) _ 2 ^ dB ext 
t nXi dt 



dJ(xj,t) 
dt 

J=r 

(7) 

where M (Xi, Xj) is the inverse of the matrix M(xi, Xj). 
Now, we can use the Euler method and the initial condi- 
tions J(xi, to) — and B ext {to) = to compute the sheet 



3 



current density J(xi,t n ) at t = t n by 

~J(Xj,t) 



J(xi,t n ) = M ~\ x i,'- 
V 3=1 



dB P 



dt 



+ J(xi,t n -i) x At. 
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where At — t n — t n _i. The same approach can be used 
for Eq|4j From the above EqjSJ we can see that the fi- 
nal sheet current density J(xi,t n ) is only dependent on 
the initial conditions and t n , not on dB ext /dt. After we 
obtain J(xi,t n ), we compute the B field outside of the 
superconducting film (z > 0) using the Biot-Savart law. 

In order to validate our numerical model, we consider 
the case of vortices induced in a superconducting strip, 
for which an analytical solution exists |22j . Applying a 
magnetic field pulse with a magnitude of B ext perpen- 
dicular to the strip surface, the magnetic flux penetrates 
the strip. We solve Eq[3] numerically and obtain the dis- 
tribution of supercurrents in the strip plotted in Figj2] 
(solid lines). The small difference in the curves arises 
from the fact that the analytical solution assumes that 
at the edge, where the magnetic flux penetrates, the sheet 
current density takes the critical value. For the numeri- 
cal computation the sheet current density is completely 
determined by Eq(3| 



0.5 



Z3 



-0.5 




— numerics 
---analytics 



0.2 



0.4 0.6 

x [units of a] 



0.8 



FIG. 2: (Color online) Distribution of the sheet current den- 
sity across the strip for B ext = — ¥ HoJc — t 0. Solid line: 
sheet current density obtained by Eqj3] Dashed line: sheet 
current density our former results [22] . 



III. Magnetic traps 

In this section, we apply Eq|4] to compute the sheet 
current density induced in the superconducting disk and 
ring by various loading fields. From the distribution of 
the sheet current density we then compute the magnetic 
field carried by the vortices by means of the Biot-Savart 
law. An intuitive picture of the vortex loading process 



can be found in 



A. Traps generated by a superconducting disk 

Consider a disk shaped superconductor in the rema- 
nent state with a magnetic loading field pulse of B ext — 
— » 1.2/irjJc - > 0. The second critical field in our calcu- 
lations is approximately I.8/in</ c . For symmetry reasons 
the radial component of the total external field is always 
zero B r (r = 0, z) = above the disk center. The z 
component of the total external field B z (r = Q,z) above 
the disk center at different heights is plotted in Fig(3] 
B z (r, z) is always zero at the surface center of the disk 
(r = 0, z = 0), because the superconductor prevents any 
magnetic field perpendicular to its surface at this point. 
The inset of FigHl shows the supercurrent distribution 
obtained from EqHl We see that the supercurrent in the 
disk is separated into two regions; one with current prop- 
agating clockwise and one with current in the counter- 
clockwise direction. Applying a small bias field on the 
order of B ext /10 perpendicular to the disk, allows us to 
cancel B z (r — 0, z = Zq) at zq and create a quadrupole 
trap for low-field seeking atoms. Fig|4] shows such a trap 
for Buas = — 0.05/io^c and zq ~ 1.2a. Increasing the 
bias field will move this trap closer to the disk surface, 
however in doing so the bias field will also induce cur- 
rents into the disk, ensuring that the field at the disk 
surface center remains 0. This puts a J c dependent limit 
on the trap depth for traps created very close to the disk 
surface. 
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FIG. 3: (Color online) Vertical component B z (r,z) of the 
magnetic field carried by vortices loaded by B ext — — > 
1.2/iqJc — > vs z above the disk center r = 0. 



The superconductor has a memory of the history of 
the external fields. Therefore it is possible to design a 
pulse sequence in which the role of the bias field is as- 
sumed by another induced supercurrent carried by the 
disk. Here we use two loading pulses with amplitudes 
of 0.8/irjJ c and — 0.4/i J c respectively to write a vortex 
pattern in the disk, which results in a self-sufficient trap 
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FIG. 4: (Color online) Equipotential lines showing a 
quadrupole trap above the disk center formed by the vortex 
field together with a perpendicular field Bb ias = — 0.05/ioJc 
The vortex loading field pulse has an amplitude of 1.2/ioJc- 
The field variation per contour line is 5 x 10~ 3 fioJ c . Inset: 
Distribution of the sheet current density along the radial di- 
rection resulting from the loading field pulse. 

above the disk center, as shown in Figj5j Choosing typi- 
cal experimental values for the critical current density of 
YBCO at liquid Nitrogen temperatures j c = IMA/cm 2 
we calculate a self-sufficient trap depth of 35/iK for a 
disk with a diameter of 1mm. This trap has the feature 
of transforming into a ring trap by applying an additional 
perpendicular bias field Bu as = 0.01/xo J c . 
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FIG. 5: (Color online) A self-sufficient trap above the disk 
center formed by the vortex field created by two loading field 
pulses with amplitudes 0.8/ioJc and — 0.4/toJc respectively. 
The field variation per contour line is 9.5 x 10~ 4 /ioJ c - 

While it is possible to split the self-sufficient trap into a 
ring, it does not necessarily imply a simple loading pro- 
cedure with a high transfer efficiency from an external 
trap via the self-sufficient trap. A more straightforward 
way of producing a ring trap from an external quadrupole 
trap is to use a first magnetic field pulse of 1.2/xp J c fol- 



lowed by a second pulse — 0.6//o</ c - Due to the nonlinear 
response of the material, this pulse sequence does not 
result in a self-sufficient trap, and instead a quadrupole 
trap can be formed far away from the chip surface by 
applying an additional bias field B^ias- By increasing 
the value of Bbi as this trap is brought closer to the chip 
surface where it eventually deforms into a ring trap. 
The perpendicular bias field lifts the field zero on the 
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FIG. 6: (Color online) Equipotential lines for an atom 
with nonzero magnetic moment. The loading field pulses are 
Bext = — > 1.2/ioJc — > — 0.6/ioJc — > 0. (a) With an ad- 
ditional bias field Bi, ias = 0.03//oJ c , two traps are obtained 
above the disk center. The field variation per contour line is 
1.5 x 10~ 3 /ioJ c - (b) Bbias = 0. 035/io J c , the two traps merge. 
The field variation per contour line is 1.54 x 10~ 3 /ioJc (c) 
Bbias = 0. 06/io Jc, a ring trap is formed. The field variation 
per contour line is 1.2 x 10 _3 /xoJ c . 

disk surface and lowers the quadrupole trap, producing 
two traps above the disk center, see Fig|6ja). Increasing 
the bias field brings the two traps closer together, until 
they merge at Bu as = 0.03/iq^c, see Figjfjfb). Increas- 
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ing Bbi as further will turn the trap into a ring trap, as 
shown in Fig{6]jc). Continuing to increase Bbias increases 
the radius of the ring trap until the trap is lowered to the 
surface. 



B. Traps generated by a superconducting ring 



In this section we consider an annular super- 
conductor. The inset of Fig|8ja) shows the su- 
percurrents carried by the vortices in the rema- 
nent state, when we apply a field pulse with a 
magnitude of 0.8/in</ c to a superconducting ring. 
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FIG. 8: (Color online) (a) A quadrupole trap formed by 
a loading field pulse B ext = — > 0.8fioJc — > and a bias 
field Bbias = — 0.028/ioJc- The field variation per contour 
line is 2 x lO~ 3 fioJc- Inset: Distribution of the sheet cur- 
rent density along the radial direction of the superconducting 
ring, (b) A ring trap formed when increasing the bias field 
to Bbias = —0.06/j,oJc- The field variation per contour line is 
5 x KTVoJo 



FIG. 7: (Color online) B z (r — 0,z) carried by the vortices 
at different z above the ring center. The loading field pulse is 
B ext = -> O.SfioJc -> 0. 



Similar to the superconducting disk, B r (r = 0,z) 
is always zero by symmetry. We plot the z compo- 
nent of the magnetic field B z (r — 0, z) generated by 
these supercurrents at different height z in FigjT] If 
an additional perpendicular field is applied to cancel 
B z (r = 0,z = zq), a quadrupole trap can be formed 
at (0, Zo). When Bbias = — 0.028^o^c a quadrupole 
trap is formed at z ~ 0.7a, as shown in Figj8ja). For 
small bias fields, increasing the magnitude of Bbias 
reduces the trap distance to the chip surface. When 
I-Bfreasl ~ 0.06/in J c , the quadrupole trap transforms 
into a ring trap, as shown in Figj8f^b). Further in- 
creasing the bias field will enlarge the radius of the 
ring trap until the trap is lowered to the surface. 
Alternatively, a ring trap may be formed by apply- 
ing two field pulses and a bias field as indicated in 
Figj9] The first pulse has a magnitude of /inJ c , the 
second has a magnitude of — 0.5[ioJ c and the bias field 
Bbias = 0.005^o^c- The ring trap of Fig[9] has a larger 
radius and distance to the surface as well as tighter 
confinement than the one shown in Fig|8[b). 
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FIG. 9: (Color online) A ring trap formed by a loading field 
sequence B ext = — > fj,o Jc — > — 0.5/io Jc — > and a bias field 
Bbias = 0.005pioJ c - The double pulse sequence increases the 
confinement of the atoms in the ring compared to the single 
pulse loading sequence. The field variation per contour line 
is 1.8 x 10" 3 (loJc- 
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IV. Conclusion 

We have investigated and designed a variety of con- 
fining geometries for ultra cold neutral atoms based on 
type-II superconducting disk and ring structures. These 
geometries include quadrupole traps, a self-sufficient trap 
and ring traps. We have shown that a quadrupole trap 
above a superconducting disk and ring can be formed 
by applying a single magnetic field loading pulse and 
adding an external bias field. Using a double pulse se- 
quence we are able to produce a self-sufficient trap above 
the disk center. We have developed a procedure for the 
formation of ring traps based on superconducting disk 
and ring structures. Furthermore we have shown an effi- 
cient loading procedure for these ring traps. Apart from 
methods outlined in this paper based on superconductors 
in the mixed state, it is possible to produce quadrupole 
traps with a ring-shaped superconductor in the Meissner 
state using flux trapped inside the ring and a bias field. 
By varying the trapped flux and bias field magnitudes 
ring traps can be produced. However, the limitations 
on the imprinted current distribution make it impossible 
to generate versatile magnetic field patterns [31] , The 
key point of using type-II superconducting structures is 
that a single microstructured element is capable of gen- 
crating tailored complex trap geometries for a variety of 
experimental implementations. This allows us to control 
the geometry of confining potentials and their parameters 
without having to fabricate a new atom chip. In addi- 
tion to this flexibility, superconductors have the intrinsic 



advantage of a low trap loss rate due to the reduced spin 
flips of the atoms close to the chip surface. The poten- 
tial contribution to the spin flip rate by vortex motion in 
the superconductor has yet to be experimentally investi- 
gated. 

An additional control over the trap geometries, com- 
plementary to the magnetic field pulses, can be obtained 
by dressing the atomic states with a radio- frequency field 
[32| [33] . This can help to lift the magnetic field zero in- 
trinsically present in quadrupole traps, thus reducing the 
Majorana spin flip rate. The low field noise in the pro- 
posed ring traps may make them an attractive candidate 
for the construction of a Sagnac-type interferometer |34) . 
A persistent flow of ultracold atoms, even an analog of the 
superconducting quantum interference device (SQUID), 
may be realized in these ring traps |35[ 136] . Taking 
into account all the intrinsic advantages discussed high- 
lights the potential of magnetic traps for neutral atoms 
based on patterned vortex distributions in superconduct- 
ing disks and rings. 
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